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ABSTRACT 



In this paper, we describe a reliable symbolic computational algorithm for inverting general cyclic 
heptadiagonal matrices by using parallel computing along with recursion. The algorithm is implementable to the 
Computer Algebra System(CAS) such as MAPLE, MATLAB and MATHEMATICA. An example is presented 
for the sake of illustration. 
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1. INTRODUCTION 

The n x n general periodic heptadiagonal matrices are takes the form: 







A 
















_ ^1 


b 2 


d 2 


a 2 


A 2 


c 2 













B 3 


b 3 


d 3 


a 3 


A 3 


c 3 











D 3 


B 3 




d 3 


a 3 








































D , 

n-3 




b n -s 




"n-3 

























-2 ^n- 


4,-1 











D nl 








-1 «„- 




A 


n 












5„ 







Where n >8. 



The inverses of cyclic heptadiagonal matrices are usually required in science and 
engineering applications, for more details (see special cases, [1-9]). The motivation of the 
current paper is to establish efficient algorithms for inverting cyclic heptadiagonal matrices of 
the form (1.1) and for solving linear systems of the form: 

Hx=r (1.2) 
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where x = (x l ,x 2 ,...,x n f ,r =(r l ,r 2 ,...,r n f . 

To the best of our knowledge, the inversion of a general cyclic heptadiagonal matrix of 
the form (1.1) has not been considered. Very recently in [9], the inversion of a general cyclic 
pentadiagonal matrix using recursion is studied without imposing any restrictive conditions 
on the elements of the matrix. Also, in this paper we are going to compute the inverse of a 
general cyclic heptadiagonal matrix of the form (1.1) without imposing any restrictive 
conditions on the elements of the matrix H in (1.1). Our approach is mainly based on getting 
the elements of the last four columns of H" 1 in suitable forms via the Doolittle LU 
factorization [10] along with parallel computation [6]. Then the elements of the remaining (n 
- 4) columns of FT 1 may be obtained using relevant recursive relations. The inversion 
algorithm of this paper is a natural generalization of the algorithm presented in [9]. The 
development of a symbolic algorithm is considered in order to remove all cases where the 
numerical algorithm fails. 

The paper is organized as follows. In Section 2, new symbolic computational algorithms, 
that will not break, is constructed. In Section 3, an illustrative example is given. Conclusions 
of the work are given in Section 4. 



2. Main results 

In this section we shall focus on the construction of new symbolic computational 
algorithms for computing the determinant and the inverse of general cyclic heptadiagonal 
matrices. The solution of cyclic heptadiagonal linear systems of the form (1.2) will be taken 
into account. Firstly we begin with computing the LU factorization of the matrix H. It is as in 
the following: 

H=LU (2.1) 

where 
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(2.2) 
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The elements in the matrices L and U in (2.2) and (2.3) satisfy: 
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We also have: 



detH = Y[ a r 
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(2.12) 
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At this point it is convenient to formulate our first result. It is a symbolic algorithm for computing 
the determinant of a cyclic heptadiagonal matrix H of the form (1.1) and can be considered as 
natural generalization of the symbolic algorithm DETCPENTA in [9]. 

Algorithm 2.1. To compute det H for the cyclic heptadiagonal matrix H in (1.1), we may proceed as 
follows: 

Step 1: Set cti=di. If ai=0 then a x = t end if. Set gi=ai, Zi=Ai, ki=A n _i/ a 1( Vi=b 1; w a =Bl, hi=a n / ai, 
Wi=Bi , f2=b 2 / c<i, e 3 =B 3 / ai, a 2 =d 2 -f2 *gi- If a 2 =0 then a 2 = t end if. Set K 2 =-ki*gi/ a 2 , v 2 =B 2 - 
f 2 *Vi, w 2 =-f 2 *Wi, h 2 =(A n -hi*gi)/ a 2/ g 2 =a 2 -f 2 *z 1/ f3=(b 3 -e 3 *gi)/a 2 , a3=d 3 -e 3 *z 1 -f3*g 2 . If a 3 =0 then 
a 3 = t end if. Set k 3 =-(ki*Zi+k 2 *g 2 )/a 3 , h 3 =-(hiZi+h 2 g 2 )/a 3 , v 3 =-e 3 *Vi-f 3 *v 2 , w 3 =-f 3 w 2 -e 3 W!. 
Step 2: Compute and simplify: 
For i from 4 to n-2 do 

e i =(B i -D i * gi. 3 / cti-s)/ ai-2 
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If a,=0 then a,= t end if 
End do 

Step 3: Compute and simplify: 
For i from 4 to n-5 do 

K i =-(k i . 3 *C i . 3 +k i . 2 *z i . 2 +k i . 1 *g i . 1 )/ di 
Wi=-(Di*Wi-3/ ai.3+ei*Wi. 2 +fi*Wi.i) 
End do 



Step 4: Compute and simplify: 
For i from 4 to n-4 do 

hi=-(h i _3*Ci_3+hi_ 2 *z i _ 2 +h i . 1 *g i _ 1 )/ a; 
Vi=-(Di*Vi-3/ a i _3+e i *v i _2+fi*v i -i) 
End do 
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Step 6: Compute det H 



The symbolic Algorithm 2.1 will be referred to as DETCHEPTA. The new algorithm DETCHEPTA 
is very useful to check the nonsingularity of the matrix H when we consider, for example, the 
solution of the cyclic heptadiagonal linear systems of the form (1.2). 

Now, when the matrix H is nonsingular, its inversion is computed as follows: 

Let H 1 =(Sy) HiJ ^ n =(Col 1 ,Col 2 ,...,Col n ) (2.14) 

Where Col m denotes the m th column of H ~ , m = 1,2,. 

Since the Doolittle LU factorization of the matrix H in (1.1) is always possible then we can use 
parallel computations to get the elements of the last five columns Co/, = (S li ,S 2J ,...,S ni ), 

i =n,n -l,n -2,n -3and« -AoiH~ l as follows [9]: 
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The remaining (n - 5) columns are obtained by using the fact H ~ l H = I n , where I n is the 
identity matrix. They are as in the following: 

Col n _ 5 =—{E n _ 2 -A n _fol n _ 4 -a n _foI n _ 3 -d n _fol n _ 2 -b n _fol nA -B n Col n ), 

(2.39) 

Col j =-^( E J+ 3 -A J+ fol J+1 ~a j+2 Col j+2 -d j+ fol j+3 -b j+ fol j+4 -B j+ fol j+5 -D j+6 Col j+b ). 

j=n-6,n-l,...,\ (2-40) 
where Ej=l if row=i. 

n-5 

Remark 2.1. Eqs. (2.39) and (2.40) suggest an additional assumption ]^[C ; ^0,, which is 

i=i 

only formal and can be omitted by introducing auxiliary parameter t in Algorithm 2.2 given 
below. 

Now we formulate a second result. It is a symbolic computational algorithm to compute 
the inverse of a general cyclic heptadiagonal matrix of the form (1.1) when it exists. 

Algorithm 2.2. To find the n x n inverse matrix of the general cyclic heptadiagonal matrix H 
in (1.1) by using the relations (2.15)_(2.40). 

INPUT: Order of the matrix n and the components, Di, Bi, b f , dj, a„ A i; Q, i=l,2, ...,n, 

where Di= D 2 =D 3 =C n -2 =C n -i = C n =0. 
OUTPUT: Inverse matrix, H 1 =(S ij \< i Jin . 

Step 1: If Q=0 for any i = 1, 2, n-5 set Q= t (t is just a symbolic name). 
Step 2: If Bi=0 for any i = 6, 7, n, B|=t. 

Step 3: Use the DETCHEPTA algorithm to check the nonsingularity of the matrix H. If the 

matrix H is singular then OUTPUT ('The matrix H is singular'); Stop. 
Step 4: For i = 1, 2, n, compute and simplify the components Sj, n , S i n _i, S i)n _ 2 , S i)n _ 3 and Si, n . 4 

of the columns Cj, j=n, n-1, n-2, n-3 and n-4, respectively, by using (2.15)-(2.38). 
Step 5: For i = 1, 2, n, compute and simplify the components Sj, n -s by using (2.39). 
Step 6: For j = n-6, n-7, 1, do 
For i = 1, 2, n, do 

Compute and simplify the components Sy by using (2.40). 



End do 
End do. 

Step 6: Substitute the actual value t = in all expressions to obtain the elements, S liS , i, j = 1, 
2, n. 

The symbolic Algorithm 2.2 will be referred to as CHINV algorithm. The Algorithms 2.2, 
2.3 and 2.2 in [9], [11] and [12], respectively, are now special cases of the CHINV algorithm. 

3. An illustrative example 

In this section we give an example for the sake of illustration. 

Example 3.1. Consider the 10 X 10 cyclic heptadiagonal linear system 
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By using the coefficient matrix of the system (3.1) and applying the CHINV algorithm, we get: 



H- 1 = 
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By simple calculation we obtain the solution of cyclic heptadiagonal linear system (3.1) 
X=[xi, x 2 , x 3 , x 4 , x 5 , x 6 , x 7 , x 8 , x 9 , x 10 ]=[l,2,3,4,5,6,7,8,9,10]. 

4. Conclusions 

In this work new recursive computational algorithms have been developed for computing the 
determinant and inverse of general cyclic heptadiagonal matrices and for solving linear systems of 
cyclic heptadiagonal type. The algorithms are reliable, computationally efficient and will not fail. The 
algorithms are natural generalizations of some algorithms in current use. 
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